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ABSTRACT 

We investigate the circumstellar dust properties of the oxygen-rich bipolar proto- 
planetary nebula IRAS 18276-1431 by means of two-dimensional radiative transfer 
simulations of the circumstellar dust shell. The model geometry is assumed to have a 
torus and an envelope which consists of a pair of bipolar lobes and a spherical AGB 
shell. The parameters of the dust and the dust shell are constrained by comparing the 
spectral energy distribution (SED) and near-infrared intensity and polarisation data 
with the models. The polarisation in the envelope reaches 50 - 60 % and is nearly 
constant in the H and K$ bands in the observations. This weak wavelength dependence 
of the polarisation can be reproduced with a grain size distribution function for the 
torus: 0.05 /zm ^ a with n (a) oc a~( p=5 - 5 ' exp (— a/a c = 0.3 /mi). The power index p 
is significantly steeper than that for interstellar dust (p ~ 3). Similar results have also 
been found in some other PPNs and suggest that mechanisms that grind down large 
particles, such as sputtering, may also have acted when the dust particles formed. The 
spectral opacity index ft is found to be 0.6±0.5 from the 760 /mi to 2.6 mm fluxes, 
which is characterised by the dust in the torus. This low value (<2) indicates the 
presence of large dust grains in the torus. We discuss two possible dust models for 
the torus. One has a size distribution function of 1.0 (im ^ a ^ a max = 5 000.0 /im 
with n{a) oc a~(p= 2 - 5 ) and the other is 1.0 /tm ^ a ^ a max = 10 000.0 /mi with 
n{a) oc a~( p=3 - 5 ) . The former has (3 of 0.633, but we are not able to find reasonable 
geometry parameters to fit the SED in the infrared. The latter has (3 of 1.12, but 
reproduces the SED better over a wide wavelength range. With this dust model, the 
geometric parameters are estimated as follows: the inner and outer radii are 30 AU 
and 1000 AU and the torus mass is 3.0 Mq. Given that the torii are generally not 
found to be rotating, a large fraction of the torus material is likely to be expanding. 
Assuming an expansion velocity of 15 kms -1 , the torus formation time and mass- loss 
rate are found to be ^300 yrs and ~T0~ 2 AfQyr -1 respectively. 

Key words: Stars: AGB and post- AGB - circumstellar matter - radiative transfer 
- polarization - individual (IRAS 18276-1431) 



1 INTRODUCTION 

IRAS 18276-1431 (hereafter 118276), also known as OH 
17.7-2.0, is an oxygen-rich proto(pre)-planetary nebula 
(PPN). This object has been well studied in maser obser- 
vations over the past three d ecades since it s discovery as a 
strong OH 1612 M Hz source (lBowersjll978l ). No SiO maser 
has been detected |Nvman et al.ll 199S ) . The intensity of the 
H2O maser dropped by a factor of a hundred during the 
period between 1985 and 1990 and is now undetectable, in- 
dicating that the ma ss-loss rate has rapidly decreased to 
below 10~ 7 M yr _1 (|Engelsll2002l ). These maser properties 



suggest that the central star evolved off the asymptotic giant 
branch (AGB) phase very recently. 

The dust she ll of this object has also been t he subject 
of several studies. ISanchez Contreras et all (|2007l . hereafter 
SC07) presented optical HST/WFPC2 images and infrared 
Kp, Lp, and Ms images from the Keck II telescope. The im- 
ages show a striking bipolar appearance with an extension of 
~2 arcsecx3 arcsec. While the extension and the appearance 
look similar in these shorter wavelength ranges, the mid- 
infrared images at 8.59 /im, 11.85 /mi, and 12.81 /im show 
an elongated shape along the polar direction with a flux 
peak towards the central star (jLagadec et al.l [201lf ) . SC07 
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also performed a one-dimensional (ID) double shell model 
that mainly fits the mid-infrared flux. They found that the 
observed fluxes between 800 /jm and 2.6 mm are roughly 
proportional to v 1 , where v is the frequency, indicating the 
opacity coefficient Q v ~ 1. Hence, they concluded that the 
dust grains with radii a > 400 /im could explain the flux 
slope at these wavelengths. 

The dust growth in the circumstellar environment is an 
interesting issue. If the ejected material radially expands at 
typical expansion velocities of ~15 kms , the dust particles 
may only grow a little and retain the properties that were set 
in the dust formation region. That is to say, particles with 
sizes up to a few micron sho uld exist. This is con firmed in 
several previous papers (e.g. I Jura fc Morrill 19851 ). Bipolar 
objects such as 118276 are expected to have an optically 
thick, long-lived dust torus or disc around the central star, 
which is t hought to hav e formed because of binary interac- 
tions (e.g. lMorris|[l987h . In such a thick dusty region, the 
dust grows in size by mutual collisions. Although the forma- 
tion of the disc/torus and the dust processes therein have 
been studied in various object classes, they are still under 
debate and present some intriguing problems. 

In the present work, we investigate the dust proper- 
ties in I18276's circumstellar dust shell (CDS) by means of 
two-dimensional (2D) radiative transfer calculations. An ap- 
proximate solution was obtained by SC07's ID SED model. 
However, since 118276 exhibits a striking bipolar lobe, we 
explore 2D modeling. Furthermore, the spectral slope of the 
(sub) millimeter flux is often used to detect large grains of 
a few hundred micron or larger. In our modeling, we en- 
counter difficulties explaining the millimeter flux excess and 
the fluxes at far-infrared or shorter wavelengths simultane- 
ously with a reasonable mass of the material in the equato- 
rial region. We will revisit the spectral slope issue. 



2 DUST SHELL MODELING 

We performed 2D radiative transfer modeling of I18276's 
CDS. Our aim is to investigate the dust properties in the 
envelope and the equatorial region. An analysis of the arc- 
like structure and the "search light" feature, which are de- 
tected in previous optical and NIR images (SC07), and the 
formation of bipolarity are beyond the scope of this work. A 
comprehensive discussion on the morphology can be found 
in SC07. Our modeling approach is bas i cally similar to ou r 
previous work ( Muraka wa et al.l l2010al ; [Murakawa 2 010bh . 
We used NIR polarimetric data obtained with VLT/NACO, 
these data are particularly useful to estimate the dust sizes. 
This is complemented by SED data collected from various 
sources. The result s of the HST and Ke ck images (SC07) and 
the VISIR image (|Lagadec et al.ll201ll ) are weighted less to 
constrain the parameters in our 2D modeling. In the follow- 
ing subsections, the archival data, model assumptions, and 
the results for the selected model are presented. 



2.1 Observational constraints 

Figure[TJ shows the SED of 118276, it has a single flux peak 
at ~25 /im and a weak 10 /im absorption feature. In this 
plot, the interstellar reddening, which is calculated to be 
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Figure 1. The observed SED of 118276, with the model results 
over plotted. The solid, dotted, and dashed-dotted curves denote 
the total flux, the scattered light, and the thermal emission, re- 
spectively. The black lines represent the model with the torus and 
envelope. The red line denotes the model with only the torus. The 
sources for the phot ometry are in or der of wavelength, from the 
optical to millimete r:lLe Be rtre ( 1987), Catalog of Infrared Obser- 
vations, Edition 5 l lGezari et al.lll999l) . 2MASS all -sky catalogue 
of po int sources, Akari/IRC mid-IR all-sky Survey ilshihara et al.l 
l2010t) . the calibrated data (the d ata ID of 10802940) of the ISO 
Short Wavelength Spectrometer iKraemer et al.ll2002r) . MSX in- 
frared point source catalogue, IRAS point source catalogue, and 
SC07. For the data where the error bars are not seen, the uncer- 
tainties are smaller than the plot symbols. 



1.6(±0.5) mag at V (SC07), is used to correct the measured 
fluxes for wavelengths shorter than 5 /xm. 

The upper panels of Fig. [2] present the NIR images, 
which are collected from the VLT/NACO data archive. The 
natural seeing was 0.5 arcsec to 0.6 arcsec in the optical dur- 
ing the observations, which is typical for the Paranal site. 
The H and Ks band data were taken with the Wollaston 
prism to measure the polarisation. The data were reduced 
follow ing the same method as before (|Murakawa fc Izumiural 
l2012i ). Standard star data were used to calibrate the surface 
brightness and the polarisation, and also served to estimate 
the size of the point spread function (PSF). The measured 
full widths at half maximum of the beam are 0.1 arcsec, 
0.058 arcsec, and 0.063 arcsec in the J, H , and Ks bands, 
respectively. The PSFs were modeled using one-dimensional 
Moffat function fitting. The PSF size 2<r and slope /3 are 
obtained to be 0.146 arcsec and 1.17 in the J-band, 0.074 
arcsec and 1.18 in the ff-band, and 0.048 arcsec and 1.03 
in the JCs-band. The model images are convolved with the 
model PSF which allows comparison with the observed im- 
ages. The intensity images show a striking bipolar appear- 
ance and look similar in the J, H, and Ks bands, consistent 
with previous observations (e.g. SC07). For the H and Ks 
band images, the polarisation vectors are overlaid. The vec- 
tors are aligned perpendicularly to the bipolar axis. The 
polarisation reaches 50 - 60 % and does not change much 
between the H and Ks bands. There is no significant differ- 
ence between the upper, brighter and lower, fainter lobes. 

In our modeling, we aim to fit the shape of the SED in- 
cluding the weak 10 /im silicate feature, the constant bipolar 
appearance and degree of polarisation in the NIR. 
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Figure 2. Left three columns: J, H, and ifg-band images. Top and bottom panels are from the VLT/NACO archived and our model 
results, respectively. The color scale bars indicate the surface brightness in Wm" 2 pm _1 arcsec — 2 . In the H and ifg-band images, the 
polarisation vector lines are overlaid. The modelled images are convolved with a one-dimensional Moffat model PSF. 



2.2 Model assumptions 

The spectral type and the total luminosity of the cen- 
tral star are estimated to be ^K5 or earlier and ~ 1.5 x 
10 4 (D kpc/4.6) 2 , respectively (|Le Bertre et all 1 19891) . Pre- 
vious OH maser phase-lag measurements yielded a distance 
in the range 3.4 to 5.4 kpc |Bowers et alj|l98ll . ll98"3h . In our 
model, the central star is assumed to have a blackbody spec- 
trum with 7000 K, which is an intermediate value of the pre- 
vious estimates of 4000 K to 10 000 K (e.g. iLe Bertre et all 
Il989h . We adopted the distance to be 3 kpc following SC07, 
yielding a luminosity of 6380 Lq. We will start with this 
value, but will eventually use 8500 Lq to better fit the SED. 

In bipolar PPNs such as 118276, it is expected that the 
CDS has an optically thick dust torus/disc in the equato- 
rial region and an envelope, which characterises the lobe 
shape. In our separate paper on the bipolar PPN IRAS 
16342-3814, we needed an addition al optically thick, geo- 
metr ically thin disc inside the torus (jMurakawa fc Izumiural 
2012). For 118276, we tried both model geometries, but 
found that the inner disc is not necessary. The mass den- 
sity form for 118276 is assumed to be p — ptorus + 
piobc + Pages • We apply a mass density form ptorus (r, z) = 
p t (?7i?torus)~ 2 exp (— z 2 /2 (Hor) 2 ^) for the torus, where 
(r, z) are the two-dimensional cylindrical coordinates, and 
the torus region is defined as Ri n ^ Vr 2 + z 2 < -Rtorus- We 
would like to stress that the dust particle motions in the 
torii of evolved stars are not well known. For example, it 
is not clear whether they are accreting onto the star, ex- 
panding, or rotating in a Keplerian fashion, or not. Hence, 



we regard the torus just as a tapered dust structure. For 
the envelope, we assume a spherical AGB shell in the outer 
part and a pair of bipolar lobes inside the AGB shell. This 
is because the mass-loss rate rapidly increases at the very 
end of the thermal pulsing AGB phase. The lobe is then 
carved out by the fast wind when the star evolves to the 
post-AGB phase. The bipolar appearance is mainly charac- 
terised by the lobe structure. The mass density forms of the 
lobe and AGB shell are given by equations (4) and (5) in 
our 116342 paper, respectively. They are determined by the 
density coefficient, p e , the density ratios ei n and e r im, for the 
inner cavity and the rim of the lobe; the lobe shape, f), the 
lobe thickness, 7, the lobe size, R\obe, and the outer radius 
Rout- The mass density coefficients, p t and p c , are deter- 
mined from the masses of the torus and envelope (i.e. lobe 
and AGB shell) assuming a gas-to-dust mass ratio of 200. 
The important parameters are the radii of the inner 7ii n and 
outer -Rtorus boundaries of the torus, the aspect ratio of the 
torus height Ho, the masses of the torus M t0 rus and the en- 
velope M onv , 7, £i n , and e r j m are determined in our modeling. 
The other parameters such as f3, Riobe, and R ou t are set to 
be fixed values. Figure[3] shows the mass density map of the 
inner 2000x3000 AU region. 

The mineralogy of dust in the CDS is not well known. 
We therefore simplify the dust model and adopt the opti- 
cal c onstants of oxygen deficient silicate l|Ossenkopf et al.l 
Il992l) . We assume different dust sizes for the torus and the 
envelope. The dust size in the envelope can be constrained 
with polarimetri c data. We assum e a KMH-like size distri- 
bution function (|Kim et al.l[l9 94): n (a) oc a~ p exp (— a/a c ) 
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Figure 3. Mass density distribution of the model in a plane cut- 
ting through the symmetry axis. The parameters are from the 
selected model (see Table [TJ. The AGB shell is not explicitly in- 
dicated, because the mass density is lower than 2xl0~ 21 gem -3 , 
but it exists outside the bipolar lobe and torus. 



Table 1. Model parameters of I18276's circumstellar dust shell. 



parameters 


values 


comments 1 




central star 


T e ff 


7000 K 


Le Bertre et al. (1989} 


L* 


8500 L Q 


adopted 2 


D 


3.0 kpc 


SC07 


R* 


4.38 x 10 13 cm 


calculated 




torus 




Rin 


30 AU 


20 - 30 


RtOTUS 


1000 AU 


700 - 1500 


H 


1.0 


0.8 - 1.2 


A^torus 


3.0 M Q 


2.5 - 4.0 


^min 


1.0 fira 


adopted 




10 000.0 fim 


> 1000.0 


^/im 


52 


calculated 




bipolar lobe and AGB shell 


-Rlobc 


3000 AU 


adopted 3 


P 


1.0 


adopted 


1 


0.75 


0.7 - 0.8 


Qn 


0.01 


< 0.1 


e rim 


10 


10 - 15 


Rout 


90 000 AU 


adopted 




1.2 M & 


1.0 - 1.5 


a c 


0.3 (im 


0.3 - 0.4 



1 Ranges give the uncertainty of the corresponding model param- 
eters. 2 Based on comparison of the SED. 3 Based on comparison 
with the intensity image. 



with a m i n of 0.05 /jm and note that a steep power index p 
of 5.5 results in a weak wavelength dependence. Only the 
cutoff size, a c , is a free parameter. The wavelength depen- 
dency of the power index p will be discussed in Sect. 13. ll 

The dust sizes in the torus can be constrained with 
the spectral flux slope a in the (sub) millimeter wavelength 
ranges. With the flux data between 761 fim. and 2.6 mm, 
a is 2.6T0.5, which yields the spectral opacity index f3 of 
0.6T0.5 (= a — 2). In principle, this allows us to determine 
a size distribution function, however, it is not a unique so- 
lution. We encounter difficulties in fitting the SED at the 
shorter wavelengths with the derived dust models. The de- 
tails of this problem will be discussed in sect . 13. 21 Instead 
of applying a dust model that can be determined using ft, 
we attempted to find one that fits the SED better over a 
wider wavelength range. Assuming an MRN-like size dis- 
tribution function (jMathis et al.lll977f l: n(a) oc a" 3 ' 5 with 
aWn ^ a < a max , a max is determined by the SED fitting. 
For a m i n , we adopt 1.0 fj,m. However, since we do not have 
enough clues to determine the value precisely, we exclude 
the a m i n value from the discussion of the dust properties in 
the torus in Sect. 13.21 

2.3 Results and comparison with the observations 

The aforementioned parameters are estimated by means 
of radiative transfer c alculations using our STSH code 
(I Murakawa et al.l l2008al ) . We first tried several parameter 
sets to get approximate solutions and to understand the pa- 
rameter dependence on the model results. Then, a few thou- 
sand parameter sets were examined by SED fitting. From 
this result, a dozen good parameter sets are chosen and 
their images are evaluated. We finally chose one that ex- 
plains the observations fairly well. We would like to stress 
that it is very difficult to find a parameter set that repro- 
duces all aspects of the observations. Because our aim is to 
investigate the dust properties, our modeling is mostly con- 
cerned with explaining the NIR polarisation and the SED 
over a wide wavelength range. The images in the optical and 
MIR are not accurately reproduced. The derived parameters 
are: i? in = 30 AU, i? t or US = 1000 AU, Af torus = 3.0 M & , 
aWn = 1.0 /im, a max = 10 000.0 /im, the optical depth in 
the torus midplane is 52 at 5.0 /jm, M cnv = 1.2 Mq, and 
a c — 0.3 fim. The viewing angle measured from the polar 
axis is 80°. Table [T] lists these parameters and uncertainties. 
We briefly summarise the selected model results in this sec- 
tion and discuss the effects of dust sizes on the model results 
in detail in sect. [3] 

The SED of the selected model is plotted with black 
lines in Fig.[T] The aperture sizes used for the photometric 
and spectroscopic data are different for the various obser- 
vations; they are 4 arcsec in radius (for 2MASS) or larger. 
Because the apparent size of the nebulosity is small and most 
flux is concentrated within ~1 arcsec, the effect of aperture 
size on the observed SED is negligible in 118276. In fact, 
the model flux derived applying a 4 arcsec aperture differs 
less than 1 % from that applying an infinite aperture from 
optical to millimetre wavelengths. The model fits the ob- 
served SED well apart from the optical. The weak 10 /im 
feature is due to a filling in the underlying absorption with 
strong thermal emission. In fact, the MIR images sh ow a 
flux peak towards the central star (La gadec et al.ll201ll ) and 
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the calculated optical depth becomes high as 41 at 10 /im. 
The red line denotes the result of the selected model but 
with only the torus, e.g. M env is simply set to be 0. It ap- 
proximately indicates the contribution from the torus in the 
selected model. In the wavelengths shorter than ~9 /im and 
between ~9 /im and ~200 /jm, the scattered light and ther- 
mal emission from the envelope, respectively, dominate. The 
SED of the selected model fits the observations comparably 
well in the NIR and longer wavelengths. However, a large 
discrepancy is seen in the optical. In some models with other 
parameter sets, e.g. aw n = 5 fim dust for the torus and a 
shallower power index of the mass density distribution in 
the envelope, we find that these models can fit the SED in 
the optical better than the selected model. However, these 
models provide poorer fits to the images or the MIR fluxes. 

The bottom panels in Fig. [2] show the model results. 
The pair of bipolar lobes and a nearly constant appearance 
from J to Ks bands are reproduced. The optical and MIR 
images are also reproduced although they are not presented. 
The bipolarity is seen in the optical, Lp, and Ms bands. On 
the other hand, the 11.85 /im image shows a single peak. 
These results ar e qualitatively consi stent with the observa- 
tions (SC07 and lLagadec et al]|201ll ). However the details of 
the surface brightness structures of the images are not accu- 
rately reproduced because the images at these wavelengths 
are not used to strictly constrain the model parameters. 

The polarisation vectors are aligned perpendicular to 
the bipolar axis in the entire nebulosity. Since a centro- 
symmetric pattern is seen in the raw images, this alignment 
is due to the PSF s moothing effect, as is also dem onstrated 
in our 116342 paper l|Murakawa fc Izumiuralfeo^ . The po- 
larisation in the bipolar lobes reaches 50 - 65 %, which is 
slightly higher than the observations. If a c — 0.4 /im is used, 
the polarisation becomes 45 - 60 %. We confirm that our 
estimation o f the par t icle si ze in the envelope is consistent 
with that bv lGledhilll (|2005h . 

We reach an important conclusion with respect to the 
particle sizes. In the envelope, the mass fraction of the 
particles with radius larger than the cutoff size (a c — 
0.3 fJ,m) is 28 %. This fraction decreases dramatically 
for larger sizes, e.g. only 1.7 % for a ^ 1.0 fim. The 
submicron sized dust dominates in the envelope. On 
the other hand, the fraction of the particles larger 
than 400 fim is high as 81 % for the torus. Hence, the 
particles in the torus are expected to be significantly 
larger than in the envelope. In Sect. l3.ll we discuss the 
effect of the power index, p, on the polarization and pro- 
vide our interpretation on the derived p of 5.5. In Sect. 13.21 
we discuss the determination of the particle sizes using the 
millimeter flux excess and the problems encountered in our 
modeling. 



3 DUST PROPERTIES 

3.1 In the bipolar lobe and the AGB shell 

In NIR polarimetric images of bipolar reflection nebulae 
such as those seen around for example evolved stars, young 
stellar objects (YSOs), and active galactic nuclei, the po- 
larisation vectors are centro-symmetrically aligned and the 
degree of linear polarisation reaches > 10 % in general. This 
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Figure 4. Plot of scattering property. The degree of polarisation 
is the maximum values of the single scattered light. Red and black 
lines denote the results of Gl and G2, respectively. 



result is explained by single scattering off submicron-sized 
dust particles. The polarisation is sensitive to the dust size, 
or to be more precise, the size parameter x = 2na/X for 
x ~ 1. In the case of YSOs, a strong wavelength dependence 
of the polarisation is seen which can be ex plained with a size 
distribution function: n (a) 



.. lHales et al.ll2006l ; 
iBeckford et~ai1l2008l : iMurakawa et al.ll2008aT ). This value for 
the exponent can be readily understood because the dust 
in the YSO envelope originates from interstellar matter. In 
contrast, weak wavelength dependencies are found in some 
evolved stars. The polarisation rang es from 35 to 60 % in 
the optical to NIR in F r osty Leo ([Dougados et alj Il990l ; 
IScarrott fc Scarrottll994h . IScarrott fc ScarrottJ ljl994h found 
a steeper power index p ranging from 5.2 to 5.5. 

We now investigate the effects of the dust sizes and 
the size distribution function on polarisation. Two mod- 
els are considered: Gl has 0.05 /im ^ a and n (a) oc 
a -5 ' 5 exp (— a/a c ), and G2 has 0.005 /im ^ a ^ a max M m 
and n (a) oc a~ 3 ' 5 . Upper limits of the dust size a c or a max 
of 0.1, 0.2, 0.3, 0.5, and 0.7 fim are examined. Figure [4] 
presents the resulting polarisation. The degree of polarisa- 
tion P max is calculated using the ratio of — S12 to Sn at 
the scattering angle that maximises the polarisation, where 
Su and 1S1 2 are the scattering matr ix elements for spher- 
ical grains |Bohren fc Huffman! H983T I. These values should 
be the upper limit of the polarisation expected in the enve- 
lope because the scattered light at various scattering angles 
is integrated in the observed data. In fact, the polarisation 
values in our model are 50 - 65 % in the H and Ks bands, 
but the above calculations give 70 - 80 %. 

G2 has a strong wavelength dependence in the opti- 
cal to the NIR, as is observed in YSO envelopes. On the 
other hand, Gl has a weak dependence in wavelength and 
dust size. The scattering properties are determined by the 
relative fractions of small and large dust. The scattering 
matrix elements S are weighted with n(a) by < S >= 
J Sn(a)da/ Jn(a)da. In G2 the scattering property is 
dominated by the dust at the maximum sizes. The polar- 
isation varies rapidly around x ~ 1, e.g. A ~ 1 /im for 



© 2002 RAS, MNRAS 000, [L>?? 



6 K. Murakawa et al. 



fflmax = 0.3 /mi. It is of interest that in fact the polarisa- 
tions cannot be simultaneously fit in the H and Ks bands 
if the G2 dust is applied in the envelope in our modeling. 
For example, whereas the a max = 0.4 /im model reproduces 
the H band polarization with a value of Ph ~ 55 %, the 
model returns a too high K band polarization of ~ 70 %. 
The 0.5 fim dust models, instead, reproduce the K band 
polarization and return Pk ~ 60 %, however they result 
in too low Ph ~ 27 %. We see large discrepancies despite a 
small difference of ~0.1 /im in the a max values. On the other 
hand, the Gl model has a larger fraction of small particles 
(x < 1). Thus, the polarisations are higher even at shorter 
wavelengths, i.e. A < 1 /mi. In contrast, lower polarisations 
are obtained at longer wavelengths (A > 1 (im) because of 
the large particles, even though the fractions are low. 

It is also useful to note that the polarisations of the Gl 
models are lower than those of G2 at A > 3 /im even with 
the same upper limit sizes. This is because dust models with 
the KMH-like size distribution function have particles with 
sizes exceeding a c . If a steeper power index, e.g.p = 5.5 > 
3.5, is assumed in the MRN-type dust model, the scattering 
properties at longer wavelengths (A > 27ra max ) consequently 
becomes similar to those with a shallower power index. We 
reach two conclusions. First, in objects that show a weak 
wavelength dependence in their envelopes such as 118276 
and Frosty Leo, a steep power index p(> 3.5) is expected in 
the size distribution function. Second, if the polarisation is 
measured over a wider wavelength range such as from the 
optical to 3 to 5 /im, which is still in the scattering regime in 
cool dust regions, the size distribution functions are better 
constrained, e.g. the lower limit and whether the MRN-type 
or the KMH-type is more favoured. 

This brings us to the more general question whether 
the dust properties depend on the objects under consid- 
eration. A NIR polarimetric survey of a dozen of PPNs 
was presented jGledhiU et al.l 120011 ). Four objects (IRAS 
17436+5003, IRAS 19114+0002, IRAS 18095+2704, and 
IRAS 19500-1709) out of seven for which the nebulosity 
is spatially resolved, show nearly constant polarisations in 
the J and K bands. Including three other objects (Frosty 
Leo, 116342, and 118276), a weak wavelength dependence on 
polarisation is detected. We thus find similar polarisation 
behaviour for different morphologies, i.e. spherical, elliptic, 
and bipolar, while all objects except for IRAS 19500-1709 
are oxygen-rich. Differences in morphology are thought to re- 
sult from the different initial stel lar mass, i.e. the higher th e 
mass, the more asymmetry fe.g. ICorradi fe Schwarzlll995| y 
These observational results suggest that the initial stellar 
mass and the dust shell morphology do not play impor- 
tant roles in determining the dust size distribution func- 
tion. The po wer index of ^3 th at is found in the interstel- 
lar medium (Ma this et al.lfl977h is thought to be the result 
of dust formation proces s of mass-losing evolved stars (e.g. 



Table 2. Dust model parameters for the torus 



or dust formation process 
iBierman fc Haiwrdll980l) 



Dominik. Gail, fc Sedlmavrl 0989) analysed the dust 



size spectrum in the dust-driven winds of carbon rich stars. 
In their models, dust particles form via nucleation at < 
2.0 R+, the size distribution function rapidly approaches 
the final distribution within ~20 R+. They found that the 
size distribution functions are proportional to a~ 5 in the 
case without drift between the gas and dust motions and to 
a~ 7 in the case with drift. Their results show that large 



ID 


Amax / P 


K 


Cext 


A^torus 




/im/ 


9 — 1 

cm g 


9 — 1 

cm g 


M 


A 


10 000.0/3.5 


1.12 1.44 


103 


3.0 2.11 


B 


5000.0/2.5 


0.633 1.18 


4.54 


68.1 2.59 



The spectral opacity index j3 and the mass absorptive opacity 
k are values at A = 1.1 mm. The value of the extinction cross 
section C GX t is for A = 5.0 fim. For the torus mass Mtorus, the 
values are obtained to keep a constant optical depth at 5.0 (im 
t — 52 (left column), i.e. C ex tM t orus is constant, and to fit the 
flux at A = 1.1 mm (right column). 

drift velocities steepen the size distribution functi on be- 
cause they can break up large pa r ticles by sputtering. IWoitkd 
120061 ) and iHofner fc Andersen! |2007h pointed out that in 
oxygen-rich winds, the radiation pressure on the dust is too 
weak to drive the wind because of low dust opacities. In 
fact, survey observations of a large number of AGB stars 
and mass-loss models have show n that expansion veloci- 
ties are lower in oxy gen-r i ch stars jGroenewegen et al.lll999l ; 



tics arc lower m oxygen-nen stars ( Urocncwcgcn ct al. tyyy; 
ISchoier fc Olofssoni l200ll ; lOlofsson et al.l l2002h . suggesting 



lower drift velocities in oxygen-rich stars than the carbon- 
rich objects. At first glance, this is in the opposite sense to 
the result that the majority of the objects where steeper 
power indices are found are oxygen-rich in the aforemen- 
tioned polarisation survey. Follow up observations of vari- 
ous objects, i.e. C-/0-rich, and various stellar luminosities, 
mass-loss rates, and evolutionary stages, and further anal- 
yses of the dust formation and stellar wind are essential. 
Although we are not in a position to identify the physical 
reasons, one may expect that shattering or sputtering works 
efficiently in the circumstellar environments in objects like 
118276. 



3.2 In the torus 

Measuring the flux excess at submillimeter to millimeter 
wavelengths provides a simple and powerful method to eval- 
uate whether large dust grains (> a f ew hundred microns) 
exist in the disc/torus region o r not dBeckwith fc Sargenj 
1 199 it iMivake fc Nakagawaf[l993l ; lDraineil2006r ). In fact. nvT 
merous previous works have found du st growth particu- 
larly in discs around young stars (e .g. ICalvet et alj |2002| ; 
iTesti et al 1 l2003l ; lLommen et ai1l2007f ) . Some intensive mod- 
eling to investigate the effects of th e dust growth on th e 
SEDs and images has been done fe.g. iD'alessio et alj|200lh . 
In our modeling of 118276, the flux from the torus dominates 
in the FIR or longer wavelengths (see Fig. [T] and Sect. 12.31) . 
In this result, the dust size in the torus was estimated to 
fit the SED in a wide wavelength range. In this section, we 
discuss dust models that fit the flux slope in the submil- 
limeter and millimeter wavelength ranges by applying the 
aforementioned analysis. 

The spectral opacity index cannot be directly deter- 
mined from observations. Rigorously speaking, it should be 
derived from the (sub) millimeter flux by solving the ra- 
diative tra nsfer problem. In special c ases, such as accre- 
tion discs (|Beckwith fc Sargenj Il99ll ). approximate solu- 
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Figure 5. (top) mass absorptive opacity ft (A) and (bottom) spec- 
tral opacity indices ft at A = 1 . 1 mm as function of the maximum 
size a max . Black and red lines denote dust models with a power 
index p of 3.5 and 2.5, respectively. The minimum size a m ; n is set 
to be 1.0 /im. In the bottom panel, the dust models with a power 
index of 4.5 (gray) and 5.5 (blue) are also indicated. In the top 
panel, the plot for the a max = 1.0 )im model is identical to the 
single sized dust with a radius of 1.0 /im. 



tions are given. Here we show a simpler prescription. Let 
us assume that the flux in that wavelength regime is dom- 
inated by emission from the torus and that this emission 
can be approximated by blackbody radiation with a rep- 
resentative dust temperature, B v (Td), the flux F v is given 
by F v w MtoTusKvBv (T d ) D~ 2 (~ x , where k„, D, and £ are 
the mass absorptive opacity, the distance to the object and 
the gas-to-dust mass ratio, respectively. If the dust emission 
at this wavelength range is assumed to be in the Rayleigh- 
Jeans limit, the blackbody function can be approximated 



by a power law, B v (Td) 



In addition, if the dust opac- 



ity has a power-law dependence on frequency, k v oc u 13 , the 
flux is given by F v oc v a with a = d\nF v /d\nu = ft + 2. 
This can be obtained when the ratio of the optically thick 
to thin emissio n, A , is ne glected in the result derived by 



iBeckwith fe Sargent! (|l99lh . 

With the observed flux data of 118276 from A = 761 /mi 



to 2.6 mm, ft is found to be 0.6T0.5 where the flux measure- 
ment errors are taken into account. Next, we consider (3 in 
the dust models. Figure[5] shows the mass absorptive opacity 
k and the spectral opacity index ft at A = 1.1 mm for sev- 
eral spherical dust models, which are calculated using Mie 
scattering theory. The examined size distribution functions 
have an MRN-like form with p = 5.5, 4.5, 3.5 (G3), and 2.5 
(G4). The minimum size a m i n is set to be 1.0 /im and the 
maximum sizes that were examined range from 1.0 fim. to 
10 000 /mi. p is constant at ~ 2 for a max ^ 100.0 /im, as is 
known for IS dust. The behaviour of P does not change much 
for a m i n < 10 /im. An anomalous P > 2 is seen at 100 /im 
^ ctmax ^ 400 /im. This is due to a sudden drop in the opac- 
ity at x ~ 1, as seen in the top panel. For a max > 400 /im, 
P monotonically decreases. The behaviour of P depends 
on the power index p. The absorptive opacity is iva 2 Q (a), 
where Q (a) is the absorptive coefficient. The fraction of dust 
opacity to grains with size a or larger, f (a), is given by 
/ (a) > (<&2 - a 3 - p ) / - t&£). This should be the 

lower limit because Q depends on a and is less than unity 
at smaller a. The larger the p value, the smaller the fraction 
of large particles. In the dust models with p = 5.5 and 4.5, 
P stays ~2 even with large a max values. In the G3 type, P 
decreases down to ~ 1. Such a lower limit is reproduced i n 
astronomical silicate particles (/3 m i n ~ 0.85. iDrainel 12006). 
A very low P < 1 in an oxygen-rich environment implies a 
lower p < 3.5, as seen in G4. We find that a dust model 
with a max = 5000 /im and p = 2.5 has P — 0.633, which 
agrees with the observations. We assume the representative 
temperature of 200 K, which actually varies between 100 K 
and 700 K in the torus midplane, a distance of 3 kpc, and 
a gas-to-dust mass ratio of 200. To fit the (sub)millimeter 
flux, the torus mass should be 2.59 M^ for the observed flux 
of 156 mjy at A = 1.1 mm (|van der Veen et al.lll995h and 
the dust opacity at this wavelength of 1.18 cm _2 g _1 . 

With this dust model, we encounter the serious problem 
that the NIR flux cannot be reasonably fit. The reason is 
that the extinction cross section of this dust model is too 
low in the NIR. It would be easier to understand this using 
our selected model as reference. The optical depth of the 
torus is found to be 52 at 5.0 /im. To fit the NIR flux with a 
dust model other than that used for the selected model, the 
optical depth of the torus should be kept constant. Namely, 
C cxt Mtorus should be constant because T tor us oc C cx tM t oru S , 
where C cx t is the extinction cross section at A = 5 /im. 
Table[2]lists P, k, C cx t, and torus masses of two dust models. 
Whereas C ox t of the dust model A (the selected model shown 
in Figure[T]), is high at 103 cm 2 g~ 1 , that for model B, which 
fits P, is only 4.54 cm 2 g _1 . The torus mass for model B is 
required to be unreliably large at 68.1 Mq. Hence, the P 
value is more likely close to unity. 

In either estimations (Mtorus = 3.0 Mq or 2.59 Mq), 
the torus is massive. Even if a binary companion with mass 
of ~1 Mq is assumed to exist in I18276's stellar system, 
it is unlikely that the derived torus has a pure Keplerian 
rotating motion. Rather a large fraction of the torus ma- 
terial should be radially expanding. A hydrodynamic sim- 
ulation of a binary interacted outflow model shows a spi- 
ral structure in AG B shells even in the models that are 
classified as bipolar (|Mastrodemos fe~M orris 1999). Similar 
features are detected in s ome AGB shells, e.g. NGC 306 8 
(M auron fc Hugginsll2006l ) and R Scl l|Maercker et alj|2012l ). 
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although these are optically thin cases. In addition, up to 
now, a Kepleri an rotating motion is d etected only for the 
Red Rectangle (|Buiarrabal et ai] 120031 ) while other objects 
are found to have expanding motions for which the velocity 
is well studied from observations. Assuming an expansion 
velocity, V cxp , of 15 kms~ , which is an intermediate value 
of 12 - 17 kms -1 jBowers et al l Il983l ; iHeske et all Il990l . 
SC07), the torus is formed in ~300 years (= -Rtorus/Kxp). 
To accumulate a mass of 3.0 Mq, the mass-loss rate should 
be high as ~10 -2 Mgyr -1 . Large mass-loss rates lasting a 
short time at the very end of the AGB phase have been 
found in some other objects (e.g. iBuiarrabal et all l200ll ; 
ISanchez Contreras et alJliooH ). 

We briefly highlight the problem of dust growth in the 
torus. In some objects, for which l arge g rains are detected 
e.g. the Red Rectangle (Jura et al. 1995; van Winckcl et al 
19981 ) and AFGL 2688 (jBieging fe Nguyen 1996; Jura et al 



20001 '). a long-lived disc is also expected to be present. This is 



a quite natural expectation because the dust particles have 
a larger chance to collide and grow in size in a disc than 
in an expanding outflow, although this is still hypothesised. 
However, as mentioned above, only for the Red Rectangle a 
stable rotating disc is detected. The Red Rectangle's Kep- 
lerian rotating disc is spatially resolved by C O emission line 
radio interferometry (|Buiarrabal et af1 l2005). These authors 
estimate a disc mass of 0.04 Mq, a radius of 560 AU, and a 
central mass of ~0.9 Mq. This system is somewhat similar 
to the solar nebula besides the stellar luminosity. Hence, the 
dust growth in Red Rectangle can be explained with theories 
developed for the sol ar-nebula (e.g. IWeidenschillind Il980l ; 
iNakagawa et al.lll986T l. On the other hand, in s ome objects 
like 118276 (this work) and IRAS 22036+5306 (|Sahai et all 
2006), in which large fractions of large grains are expected 
in the torus, the torus mass is very high (> 1 Mq) and 
expansion motions are detected. Further studies, in which 
the dust growth theory should be optimised for the circum- 
stellar environment around evolved stars, are encouraged to 
better understand such massive large particles in the tori of 
these objects. 



4 CONCLUSION 

We modeled the CDS of the oxygen-rich bipolar PPN IRAS 
18276-1431 by means of two-dimensional dust radiative 
transfer calculations. The previously observed SED and NIR 
polarimetric images were used to constrain the physical pa- 
rameters. The primary aim of our work was to investigate 
the dust properties. 

The polarisation reaches 50 - 60 % in the bipolar lobes 
and is almost constant in the H and Ks bands. This weak 
wavelength dependence and high values can be reproduced 
by dust models with a steeper power index of 5.5 of the size 
distribution function than interstellar dust. If the dust in 
I18276's envelope has conserved its properties since its for- 
mation during the AGB phase, mechanisms that break up 
large grains and increase the fraction of small particles, such 
as sputtering, may have also acted. Polarization measure- 
ments at a large wavelength range for various object classes 
of C-rich and O-rich, and AGB stars to PNs will provide 
better constraints of the size distribution functions and will 



help our understanding of the dust formation and process 
in circumstellar environments. 

We also investigated the dust properties in the torus. 
Two candidate dust models were examined. One has a size 
distribution function of n (a) oc a -2,5 and a max = 5 000 /im 
and the spectral opacity index /3 of 0.633. This model agrees 
with the observed /3 of 0.6 from the 760 ^tm to 2.6 mm 
fluxes, however, it reproduces too low infrared flux be- 
cause the opacity is too small. The other dust model has 
n (a) oc a" 3 ' 5 and a max = 10 000 /im and fits the SED over a 
wide wavelength range besides the (sub)millimeter opacity 
slope (f3 = 1.12). The torus masses are large for both mod- 
els; 2.59 Mq (former) and 3.0 Mq (latter). In this massive 
torus, the material is probably radially expanding rather 
than purely rotating around the central star, as is expected 
in the previous CO emission line observations. If the expan- 
sion velocity of 15 kms -1 , which is an intermediate value 
of the previous estimations, is assumed, the torus has been 
formed in ~300 yrs with a mass loss rate of ~10 -2 MQyr -1 . 
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